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ABSTRACT 

The ultimate fates of binary companions to stars (including whether the companion 
survives and the final orbit of the binary) are of interest in light of an increasing number 
of recently discovered, low-mass companions to white dwarfs (WDs). In this Letter, we 
study the evolution of a two-body system wherein the orbit adjusts due to structural 
changes in the primary, dissipation of orbital energy via tides, and mass loss during the 
giant phases; previous studies have not incorporated changes in the primary's spin. For 
companions ranging from Jupiter's mass to ^0.3 Mq and primaries ranging from 1— 
3 Mq, we determine the minimum initial scmimajor axis required for the companion to 
avoid cngulfment by the primary during post-main-sequence evolution, and highlight 
the implications for the ultimate survival of the known exoplanets. We present regions 
in secondary mass and orbital period space where an engulfed companion might be 
expected to survive the common envelope phase (CEP), and compare with known M 
dwarf+WD short-period binaries. Finally, we note that engulfed Earth-like planets 
cannot survive a CEP. Detection of a first-generation terrestrial planet in the white 
dwarf habitable zone requires scattering from a several- AU orbit to a high-eccentricity 
orbit (with a periastron of ~i?©) from which it is damped into a circular orbit via 
tidal friction, possibly rendering it an uninhabitable, charred ember. 

Key words: stars: white dwarfs - stars: AGB and post-AGB - binaries: close - stars: 
low-mass - stars: late-type 



1 INTRODUCTION 

The ultimate fate of low-mass companions to main-sequence 
(MS) stars is of interest as substellar and stellar com- 



panions to intermediate- mass stars are plentiful (Duquen 



noy and Mayor] [1991) |Raghavan e!Ta71,|2uTol |Wright et ~ 
20111. Observationally, evidence for companions to evolved 



stars is varied. In the early subgiant phase, a number of 
giant planets and giant planet candidates have been de- 
tected (Johnson et al. 2006 20111. At the end of post- 



MS evolution, low-mass companions have been found in 
post-common-envelope, short-period orbits around subd- 



warfs and white dwarfs (see 


Maxted et al.| 2006 |Silvestri 


et al.||2007| |Charpinet et al.||2011| |Liu et al.||2012| |Rebassa- 


Mansergas et al.|2012 


also Sec ; 


J 1, and in long-period orbits 


around white dwarfs 


Farihi et al.|2005 2006 2012). 



During post-MS evolution, dynamical interactions in- 
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duced by radial expansion of the primary, strong mass loss 
via stellar winds, and tidal interactions can occur. Such pro- 
cesses have been previously studied for two-body systems 



(Carlberg et al. 2009 Villaver and Livio 2009 Nordhaus 



et al.||2010| |Mustill and Villaver||2012[) and many-body sys 



terns (|Veras et al.||2011||Veras and Tout||2012| |Kratter and 



Perets 2012 Perets and Kratter 2012). In this Letter, we 



focus on two-body interactions. We improve upon previous 
studies by including changes in the rotation rates of both 
bodies due to tidal dissipation and a varying moment of in- 
ertia. During the subgiant and giant phases, the growth of 
the moment of inertia of the primary can cause companions 
within a few AU to switch from stable orbital configurations 
to being on the unstable side of the inner co-rotation point, 
thereby facilitating a plunge into their host stars ( |Spiegel| 



2012|. On the other hand, by tracking the spin-up of the 
primary, we find that higher-mass companions can retard 
their infall — an effect not incorporated in previous studies. 



For example, although Mustill and Villaver (20121 include 



the evolution of the planet's spin (which is negligible at sep- 
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arations of a few AlQ, they neglect the changes in the pri- 
mary's spin that drive the system's tidal evolution. This is 
reasonable for the very-low-mass companions they consider, 
but calculating changes in the primary's spin is necessary for 
following the orbital evolution for higher-mass companions. 

The structure of this paper is as follows. In Section 2, 
we describe the physics and assumptions of our approach. 
In Section 3, we present the minimum semimajor axis nec- 
essary to avoid engulfment as a function of companion mass 
and zero-age-main-sequence (ZAMS) primary mass. For en- 
gulfed companions, we estimate whether or not the compan- 
ion survives the common envelope and the corresponding or- 
bital period at which it emerges. In Section 4, we comment 
on the implications of finding a first-generation, terrestrial 
planet in the white dwarf habitable zone in the context of 
binary evolution. We conclude in Section 5. 



2 TIDAL DISSIPATION AND MASS LOSS 

To study the evolution of our two-body system, we employ 
a tidal interaction model that couples the mass and radius 
of the primary star with the orbit of a secondary body. Ob- 
servational constraints on tides are difficult to achieve. IZahnl 



( 1977 1989 1 proposed a tidal theory based on turbulent vis- 



cosity, which we adopt in this work. The theory was tested 
by introducing a dimensionless parameter / that was then 
calibrated using eccentricity measurements of a sample of 



post-main-sequence binaries ( Verbunt and Phinney 1995 1 



Based on the divide between observed circularized and non- 
circularized systems, it has been argued that / is constant 
and equal to unitjj^] ( |Verbunt and Phinney|[T995 ) . 

The evolution of the semimajor axis, a, is as follows: 



da 
~dt 



da\ 

dt Aides + 



(1) 



where the first term represents the change in semimajor axis 
due to tidal dissipation (described below) and the second 
term represents the adiabatic change due to mass loss (i.e., 
(da/dt) mass -i oss w —aMt/M*). Note that, while tidal dis- 
sipation occurs in both bodies, we assume that only the 
primary loses mass. Since this material is lost from the sys- 



tem (Spiegel and Madhusudhan 20121, the second term is 



positive and acts to widen the orbit. We take the orbit to be 



circular and refer the reader to Socrates et al. (2012a I for 



intricacies in modeling tidal dissipation in eccentric systems. 

As the system moves toward synchronization, the 
change in rotation rate of each body is: 



dt 



V dt 



+ 



V dt 



(2) 



where fi* and Q c are the spin rates of the primary and com- 
panion respectively. The first term represents the change in 
the spin rate of a body due to tidal dissipation in its interior. 
The second term affects the spin rate due to changes in the 



1 A jovian companion will experience a change in spin of < 1 
part in 10 3 X (a/1 AU) 6 during a Hubble time due to the stellar 
tide raised on it, where a is the orbital separation. 

2 For a detailed description of the relation between / and the tidal 
quality factor Q'^ jGoldreich and Soter|1966|l and the implications 
of the claim that / = 1 see §2.3 and 4.1 of jNordhaus et al. I 2010 1 . 



moment of inertia (Mol). In this work, we assume that the 
companion's mass and moment of inertia do not change. For 
a spin angular momentum loss rate of Si, = (2/3)M*i? 2 f2 + 
from the primary, we have that: 



dt 



2 

3e7i 



Mi, 



(3) 



where 7* is the primary's Mol, a* = J*/(M*_Rj), and M* 
and 7* are the time rates of change of the primary's mass 
and Mol|f]7* is calculated from the stellar structure at each 
point of the evolution. Our stellar models were computed 
using the Modules for Experiments in Stellar Astrophysics 
(MESA) code ( jPaxton et aT]|2l)Tl"} and produce ~0. 55-0.9- 
M e WDs for 1-3-M© ZAMS progenitors. 



According to the viscous tide model of Zahn ( 1977 



19891, the change in semimajor axis due to dissipation of 



orbital energy in the envelope of the primary is given as 



12afc 2 ,*/ (M^\ (Mc 



T*,conv 

>' I — I f l 



Mi, 
n 



Mo 
M. 



(4) 



where n is the orbital mean motion, fe,* is the primary's 
tidal Love number which we assume to be unity. The 
convective time of the primary is taken to be r* iConv = 
(M* i envTi! 2 ./T/*) 1 ' /3 such that is the luminosity of the gi- 
ant and M ijCnv is the mass of the convective envelope. The 
change in rotation rate of the primary is then: 



dQ.i, 
~dt 



6nk 2 ,*f (M, 



M, 



n 



Mc 
Mi, 

-a, 



(5) 



The equation for Cl c can be obtained by reversing the * and 
c subscripts in Eq.[5] although at a > 1 AU, Cl c is negligible. 



3 PERIOD GAPS 

Due to the combined effects of tidally-induced orbital de- 
cay and mass-loss-induced orbital expanion, there should be 
a period gap in the distribution of low-mass companions to 
white dwarfs. A cartoon schematic of this gap is presented in 
Fig. [T] The outer edge of the gap for each binary configura- 
tion is given by the final orbital separation of the companion 
that escapes engulfment (blue circle in Fig. [I]). Since tidal 
torques drop off as a large negative power of orbital sepa- 
ration, most companions that escape engulfment experience 
essentially no tidal interactions. Therefore, the outer bound- 
ary of the gap can be approximated aS CLj crit 
where ai, cr i t is the minimum semi-major axis that escapes 
engulfment]^] M*,i is the primary's ZAMS mass and M w d is 



3 If the wind is strongly magnetically coupled to the primary, 
the Coriolis torques could render the situation somewhat more 
complicated than the model presented here. 

4 We define the companion to have been engulfed when it comes 
into Roche contact with the primary ( Kopal 1959 Eg gleton|1983"l l 
and assume that its mass remains constant during the CEP 
( Maxted et al.|2006| [Passy et al.|2012b i. 
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Remnant System 



Orbit Expands 



The Gap 



Figure 1. The period gap. A companion (blue circle) orbiting just 
exterior to the critical initial semimajor axis escapes engulfment 
such that the orbit expands. A companion (red circle) orbiting 
just interior to the critical initial semimajor axis is engulfed such 
that the orbit shrinks during the CE phase. 



the mass of the emergent white dwarf. The inner edge is de- 
termined by the fate of the engulfed companion (red circle 
in Fig. [I]) which will either emerge in a post-CE short-period 
orbit, or be destroyed during the CEP. 

The inner edge of the period gap is presented in Fig. [2] 
The "allowed region" is marked in magenta. The dinner 
boundaries (orange dash-dot lines) are determined by cal- 
culating the semimajor axis where a fraction qce ( |Livio| 



and Soker 1988 Soker 2013) of the liberated orbital en- 



ergy during inspiral is sufficient to eject the envelope as- 
suming that the envelope binding energy, 25bmd = 10 46 erg 
and the mass of the WD is Mwd = 0.5 Mq. An additional 
constraint for surviving the CE phase is that the compan- 
ion does not tidally disrupt during inspiral (Nordha us and| 



Blackman 2006 Nordhaus et al. 2011| Spiegel 20121. This 



cishred boundary (yellow dash-dot line) is shown in Fig. [5] 
Additionally, the properties of the known short-period, low- 
mass companions to white dwarfs are marked by circles and 
diamonds in Fig. [2] Attempts have been made to calculate 
qce from 3-D simulations (Rick er~and Taam||2012| 
et al.||2012a l, and from observations (Zorotovic et al 



Passy| 



2010 



De Marco et al.|2011 1 but the large range of scales involved 
in a CE-inspiral makes this a difficult calculation and there 
is, as yet, no c onsensus. In Fig. [2} we present amnor for «ce 
values of 0.25 ( |Zorotovic et al.|2010[ ) and 1. 

Figure [3] shows the minimum initial semimajor axis re- 
quired to escape engulfment as a function of primary and 
companion masses. Employing t he |Verbunt and Phinney] 
(/ = 1) calibration of the |Zahn| ( |1977| |1989| > tidal formal 
ism leads to very strong tides during the post-M S, with Q+ 
values as low as ~10 1 -10 3 at the time of plunge (Nordhaus 
|et al.|[2~010| ). The values of the contours in Fig. |3| are sen- 
sitive to the maximum radius of the stellar model; changes 
in the stellar model that influence the maximum radius cor- 
respondingly influence a^crit- The discovery of companions 
near the outer edge of the gap would constrain the nature 
of tidal dissipation and the late stages of stellar evolution. 

Note that in Fig. [3} there is a decrease in the semimajor 
axis required to escape engulfment for companions greater 
than ~60 Mj. This is because at equal separation from 
the star, a more massive companion torques the star more 
strongly in proportion to the square of its mass. As the star 
ascends the asymptotic giant branch (AGB), a sufficiently 
massive companion can halt its infall by synchronizing the 
primary. Therefore, for a very massive companion to be en- 




4 5 6 8 10 



30 40 70 100 



Figure 2. The inner edge of the period gap. Objects in the ma- 
genta region can survive CE evolution and emerge as short-period 
companions to WDs. The dinner boundaries are determined by 
assuming that qce of the liberated orbital energy during inspi- 
ral is used to unbind the CE. The a s hred boundary marks the 
location where companions tidally disrupt. The known low-mass- 
short-period companions to WDs are marked with blue circles 
jParsons et al.|2012b|c|a||Haefner et al.|2004||Parsons et al.|2010| 



et al.|2012 Pyrzas et al.|2009 


|Nebot Gomez-Moran et al.|2009 


van den Besselaar et al.|2007| 


Maxted et al.|2004| |0'Brien et al.| 



Greyscale and marker size indicate companion mass. 



gulfed, its initial separation needs to be smaller than that of 
a less massive companion. Extremely massive companions 
can tidally transfer enough angular momentum to signifi- 
cantly spin up the primary. In such cases, enhanced mass- 
loss and deformation of the primary may occur. Since these 
processes are not modeled in our calculations, we exclude 
this region of parameter space (shown in grey in Fig. [3]l . 

Figure [4] shows the fates of the known exoplanets with 
masses greater than Jupiter's and separations greater than 



1 AU (data from the Open Exoplanet Catalogue: Rein 2012 1. 
For each system, we calculate the joint stellar and tidal evo- 
lution to determine whether the planet will be engulfed in 
the post-MS stages. For those objects that avoid engulfment, 
we use the WD initial-final mass relation of Catalan et al. 



( 2008 1 to determine the post-MS orbit expansion. No known 
systems will evolve to have circum-WD planets within ~6.4 
AU. Relatively few known exoplanets will survive to orbit 
a WD. Predictions of the demographics of WD planets will 
improve with a better accounting of the distribution of MS 
planets in the 3-10 AU range, which should be possible with 
micro-lensing studies with WFIRST jGreen et al.|2012 |. 



4 ARE THERE HABITABLE EARTHS 
AROUND WHITE DWARFS? 

The annulus around a white dwarf that is amenable to habit- 
able climates ( |Kasting et al.|19 93 ; Spicg eTet al.|200"8| might 
be at an orbital separation of ~1 7?0 (Agol 2011| |Fossati| 



et al. 2012 Loeb and Maoz 2013). How could a terrestrial 



planet end up in such a potentially habitable orbit? 

There are two ways that a first-generation planet could 
end up in a close orbit around a white dwarf. It either (a) 
survives a common envelope stage with a sub-giant, red- 
giant, or AGB star, in which the companion is engulfed and 
inspirals inward until its plunge is arrested before tidal dis- 
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Figure 3. The minimum separation required to escape engulfment (ai jCr i t ) as a function of companion mass and ZAMS primary mass 
under the assumption that / = 1 (Sec. pj. The outer bound of the gap can be estimated as ai iCr it X (Af*/AfvvD)i where M t is the 
primary's ZAMS mass (the ordinate in this figure) and Mwd is the remnant white dwarf mass. This estimate suggests that the outer 
boundary of the gap ranges from ~4-40 AU for the primaries and companions considered in this figure. Regions where the companion's 
mass exceeds one tenth of the primary's ZAMS mass are shown in grey and excluded from the calculations. 




Figure 4. The fates of known planetary systems. For each of 
the ~300 known exoplanets with masses at least Jupiter's and 
orbital separations at least 1 AU, we evaluate, according to Fig. [3] 
whether the object will be engulfed by its star (red dots) or escape 
cnglufmcnt (blue dots indicate initial separation while green dots 
indicate final circum-WD separations). Current and final orbital 
separations of Jupiter are shown with magneta and yellow stars, 
respectively. Among the currently known planetary systems, none 
will evolve to have circum-WD planets within ~6.4 AU. 



ruption very close to the degenerate core (at distances of 
Ss 0.3-Rq), or (b) migrates inward after the white dwarf has 
finished forming and expelled its envelope. 

Scenario (a) — CE survival — is implausible for a com- 
panion that is less than ~6 Mj (see Fig. because such 
low-mass companions lack enough energy to unbind even 
an evolved star's envelope, a necssary but not sufficient 



2012|. Furthermore, the inspiral accelerates as the compan- 
ion moves inward ( Nordhaus and Blackman|2006 Nordhaus 



et al 



20071, which means that for the star to unbind its 
own envelope right when the companion is at 1 Rq would 
require extreme fine tuning. Even if a low-mass companion 
were to arrest its inspiral at ~1 Rq, the high temperatures 
that it would encounter deep in the envelope of the primary 
(~10 6 K) might pose a severe threat to its subsequent hab- 



itability or to its very survival (Villaver and Livio 2007). 

Scenario (b) — post- WD orbital evolution — might oc- 
cur, but could render planets that experience such a pro- 
cess uninhabitable by life as we know it. Indeed, |Zucker-| 



man et al. (20101 and others have found evidence of tidally 
shredded asteroids accreting metals onto WDs, indicating 
that processes can occur in post-main-sequence planetary 
systems that excite extreme eccentricities among low-mass 
particles. Via the Kozai mechanism (Kozai 1962; Fabrycky 



and Tremaine|2007| |Katz et al.||2011| |Socrates et al.||2012b| 
Naoz et al.||2012] |Shappee and Thompson|2012[ ) , a massive 
outer body can drive an inner companion to arbitrarily high 
eccentricities, such that asteroids tidally disrupt very near 
the WD. Similar processes certainly might drive an Earth- 
mass planet from a several-AU orbit (far enough out that 
it avoided being engulfed during the AGB phase) to a high- 
eccentricity orbit with a periastron of 0.5 Rq, from which 
tidal friction would damp it to a 1-Rq circular orbit. How- 
ever, if this happens, a large amount of orbital energy must 
be dissipated as heat: 



-3xl0 42 ergsx 



0.5M©/ \M, 



Qflnal 
R . 



(6) 



condition for CE survival (Nordhaus et al. 2010 Spiegel 



where M p is the planet's mass. Dissipation of this heat 
over a circularization timescale of t c i TC ~ 10 6 years cor- 
responds to a tidal heat flux of ~(2 x 10 7 W m -2 ) x 
(tcirc/lMyr) -1 , or an average cooling temperature of 
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-1/4 



-(4000 K) x (i c ir C /lMyr) 
for the habitability of such a planet (see Barnes et al.||2012 



which could be catastrophic 



and |Barnes fc~ Rcllcr 201J2] who considered the adverse ef- 
fects on habitability of much lower tidal power). Such a mi- 
grated terrestrial planet might exist as a charred ember in 
the habitable zone. Though late-time delivery of volatiles 
might reintroduce water to a desiccated WDHZ planet, the 
near-daily cometary impacts on the Sun ( |Ma rsdcn 2005 I im- 
ply that an Earth at a 1-_Rq orbital radius could be subject 
to cometary hits at ~(i?e / R®) 2 times this rate or once per 
several decades. These impacts would have ~100 times the 
specific energy of what wiped out the dinosaurs. If large, 
these impactors could sterilize the worlds; if small, they 
would deliver very little water. 



5 CONCLUSIONS 

We have studied the effect of post-MS evolution on the orbits 
of substellar and M-dwarf companions. For each combina- 
tion of primary and companion mass, we evolve the system 
from the ZAMS through the end of the post-MS. The evo- 
lution is governed by tidal dissipation, mass loss from the 
system and structural changes in the primary. For each bi- 
nary configuration, we determine the initial semimajor axis 
required to escape engulfment. We note that terrestrial plan- 
ets cannot survive engulfment. Therefore, a first generation 
terrestrial planet in the WDHZ must have experienced an 
enormous tidal flux and, though formally habitable, might 
be inhospitable. No currently known planets will evolve to 
have circum-WD orbits within ~6.4 AU. 
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